Measurements of strongly-anisotropic ^-factors for spins in single quantum states 
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We have measured the full angular dependence, as a function of the direction of magnetic field, 
for the Zeeman splitting of individual energy states in copper nanoparticles. The ^-factors for spin 
splitting are highly anisotropic, with angular variations as large as a factor of five. The angular 
dependence fits well to ellipsoids. Both the principal- axis directions and ^-factor magnitudes vary 
between different energy levels within one nanoparticle. The variations agree quantitatively with 
random-matrix theory predictions which incorporate spin-orbit coupling. 

PACS numbers: 73.22.-f, 72.25.Rb, 73.23.Hk 



Strategies for manipulating electron spins are being 
pursued in many fieldsH, from spin-dependent tunneling-, 
to magnetic semiconductors^, to quantum computations . 
An issue central to all of these fields is that a spin is often 
not an independent variable. Interactions with the envi- 
ronment, for instance coupling to orbital degrees of free- 
dom, can limit spin coherence times and alter responses 
to magnetic fields. Here we study the effects of spin- 
orbit coupling on spins in quantum dots by exploring in 
detail how the spin properties depend on the direction of 
an applied magnetic field. Specifically, we measure the 
^-factor for spin Zeeman splitting, the sensitivity with 
which an energy level shifts with magnetic field. The 
angular dependence of the ^-factor has not previously 
been measured for quantum dots formed from semicon- 
ductor 2-dimensional electron gases, because out-of-plane 
magnetic-field components produce strong orbital effects 
that make spin energies difficult to observe; we avoid this 
problem by studying energy levels in copper nanoparti- 
cles. We find that the g-f actors for each quantum state 
can have striking anisotropics, varying by up to a factor 
of five depending on field direction. The full dependence 
on field angle is ellipsoidal .-.as expected from the tensor 
properties of the g-factorltl In addition, the ^-tensors 
exhibit large variations between different energy eigen- 
states in the same nanoparticle, with differences both in 
the direction of maximum field sensitivity and in mag- 
nitude. The variations between levels suggest that the 
strong anisotropics originate from the intrinsic quantum- 
mechanical fluctuations of mesoscopic wavof junctions, an 
effect predicted by random-matrix theoriesQfl. 

In bulk systems and for impurities in solids, it is well 
understood that spin-orbit coupling may cause ^-factors 
to exhibit anisotropics associated with directions of the 
crystal latticeErB. However™ j-jt he results of the recent 
random-matrix-theory workQH suggest that a different 
phenomenon may enhance the effects of spin-orbit cou- 
pling in nanoscale grains. The idea can be understood 
by thinking about the wavefunction for an electron in 
a metal nanoparticle, modeled as an electron in a box 
with an irregular boundary. The standing wave corre- 
sponding to the orbital amplitude of the wavefunction 
will oscillate strongly and effectively randomly as a func- 



tion of position, and this complicated wave pattern will 
vary from quantum state to quantum state. In the pres- 
ence of spin-orbit coupling, the variations present in the 
orbital wavefunction are reflected also in the spin prop- 
erties of an electronic state; in particular, they can cause 
the ^-factor to depend on the orientation of an applied 
magnetic field (B)Q. The most general dependence on 
the field direction permitted by symmetry for a spin- 1/2 
particle allows the ^-factor for an energy level /i to be 
written in the form 



g»{iS) = {glBl+ g lBl+ g lBl)^ m . 



(1) 



where #2, and gs are the ^-factors along mutually- 
orthogonal principal-axis directions (with the convention 
9\ > 92 > #3), and £?2, and ij^are the magnetic 



field components in these direct ionstffl. Pictorially, this 
means that magnitude of the ^-factor can be drawn as 
an ellipsoid. In the absence of spin-orbit coupling, the g- 
factor for spin splitting should be isotropic. However, in 
the presence of spin-orbit coupling, random-matrix the- 
ory predicts that the three principal-axis g-factors for 
each level can differ strongly, and the principal- axis di- 
rections can vary randomly from energy level to energy 
level within the same nanoparticle. 

We have probed the spin properties of individual quan- 
tum states by using electron tunneling to measure en- 
ergy levels in single copper nanoparticles. We studied 
in detail two devices containing a nanoparticle with ap- 
proximately hemispherical shape and a diameter of ~15 
nm, based on measured mean level spacings of 0.10 meV 
(Cu#l) and 0.12 meV (Cu#2). The nanoparticle is con- 
nected to two aluminum electrodes by aluminum oxide 
tunnel junctions, as shown schematically in the inset of 
Fig. 1. The procedures for fabricating and characterising 
the devices have been described in detail previouslyBEl 
We use electron-beam lithography and reactive ion etch- 
ing to create a 5-10 nm diameter hole in a silicon nitride 
membrane. A thick Al layer is evaporated on one side 
of the sample (the bottom side in the inset to Fig. 1) 
to fill the hole, and then this is oxidized in 50 mT O2 
for three minutes to form an AI2O3 tunnel barrier. We 
form the Cu nanoparticles on top of the membrane by 
self-assembly during the evaporation of 10 A of Cu. The 




0.2 0.4 0.6 0.8 



E (meV) 

FIG. 1. Magnetic field dependence of quantum states in 
copper nanoparticle Cu#l. Main panel: A color-scale plot 
of the differential conductance (dl/dV) versus energy (E) as 
a function of magnetic field (B x ). Yellow corresponds to 2 
mS, and black to zero. Values between 2 mS and the max- 
imum conductance of 3 mS are set to white. The energy is 
equal to the value of source-drain voltage multiplied by the 
capacitance factor 12 eC\/(C\ + C2) — 0.66e. The complex 
evolution of the levels with B x is due to spin-orbit coupling 
between spin-up and spin-down quantum states. The inset 
shows a sample schematic. A small hole in a Si^Na mem- 
brane is used to make electrical contact to a single copper 
nanoparticle. AI2O3 tunnel junctions (not shown) lie between 
the nanoparticle and aluminum electrodes. 

tunnel barrier on top of the nanoparticle is then made by 
depositing 11 A of AI2O3 by electron-beam evaporation. 
Finally, we deposit a thick layer of Al to form the top 
electrode. We measure current- volt age (I-V) curves for 
electrons tunneling from one electrode to the other via 
the nanoparticle. The measurements are conducted using 
a dilution refrigerator with an electron base temperature 
of approximately 40 mK. Once V exceeds a_ threshold 
associated with the Coulomb-blockade effecttil, tunnel- 
ing via the discrete quantum states in the nanoparticle 
causes the I-V curve to take the form of a sequence of 
small steps, allowing a direct measurement of the spec- 
trum of "electron-in-a-box" statesE3. 

In Fig. 1 we plot, using a color-scale, differential con- 
ductance (dl/dV) versus energy (E) curves for differ- 
ent values of magnetic field in the x-direction (defined 
in the Fig. 1 inset), showing how the energy levels in 
Cu#l evolve with B x . In this figure, dl/dV is obtained 
by numerical differentiating the I-V curves, and the en- 
ergy is determined by scaling the source-drain voltage to 
correct for the capacitive division jof voltage across the 
two tunnel junction in the deviceLL3. As has been ob- 
served previouslyli3li3, the energy levels exhibit two- fold 
Kramers' degeneracy at B=0 associated with the electron 
spin, and they undergo Zeeman splitting as a function of 
increasing B. We determine the ^-factor from the magni- 



tude of the energy splitting between Zeeman-split states, 
in the small- linear regime: #=Ai£(B)/(/ZB|B|). The 
data in Fig. 1 exhibit different values of the ^-factor for 
different quantum states and a rich variety of avoided 
level crossings as a function of increasing B. The role 
of spin-orbit coupling in producing these two effects has 
been discussed previouslyLJ'Hllj. 

Our focus in this letter is the unexplored question of 
how the ^-factors for the individual quantum-dot states 
depend on the direction of the applied magnetic field. 
We control the field magnitude and direction using a 3- 
coil magnet which can produce 1 T in any direction, or 
up to 7 T along one axistJ. We first verified that the 
tunneling spectra were identical for B and — B, as ex- 
pected due to time-reversal symmetry. Thereafter, we 
took detailed measurements only for B x >0. For each 
sample, we measured the Zeeman splitting as a function 
of |B| along_62 different directions selected to span this 
hemispheretj. The values of the ^-factors showed strong 
anisotropics, described well by an ellipsoid for each en- 
ergy level, as demonstrated below. After fitting for the 
directions of the three principal axes of the g-tensor (hi, 
fi2, hs) for each state, we took detailed field scans along 
these directions to determine precisely the principal-axis 
^-factors, and along arcs of constant field to carefully 
measure the variation of g with field direction. The long 
measurement time (1-2 weeks for each sample) requires 
exceedingly good sample stability. 

The directional dependence of g for the first quantum 
state in sample Cu#l is shown in Fig. 2. In Fig. 2(a), 
dl/dV vs. E is plotted as a function of the angle of 
the magnetic field (with magnitude 700 mT), rotating 
the field vector around the hs principal axis for the first 
quantum state from an initial orientation along J12, where 
g = 0.51±0.05, to a final orientation along hi, where g 
= 0.86±0.08. The fitted ^-factors are in good agreement 

with the form g(a) = \J g\sm 2 (a) + g2Cos 2 (a) expected 
for an ellipsoid. The data corresponding to a rotation 
from hs (where g = 0.38±0.05) to hi about the hi axis 
are shown in plots (c-d), and for rotation from hs to 
ft 2 i n ( e "f)- The ^-factors along the hi and hs axes for 
this quantum state differ by a factor of 2.3. Fig. 2(a) 
also displays the existence of energy-level to energy-level 
variations in spin properties. As a increases from to 
7r/2, the splitting of the first quantum state grows while 
the splitting of the fourth Zeeman-pair decreases. The 
other electronic states in Fig. 2 show more complicated 
behavior due to the alignment of their principal axes in 
other directions in space. We obtain good fits to ellip- 
soids for all of the measured states, with the exception 
of those that undergo an avoided crossing with a neigh- 
boring state at a very low value of field, so that the 
^-factor cannot be determined accurately (states 6 and 
7 of Cu#l and states I, 6, and 7 of Cu#2). The fitted 
parameters of the measured ^-factor ellipsoids are listed 
in Table 1. There are large variations between different 
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FIG. 2. Dependence of Zeeman splitting on magnetic- field 
direction in sample Cu#l. (a), (c), (e) Differential conduc- 
tance (dl/dV) as a function of energy (E) and the angle of the 
magnetic field (a, /3, or 7, respectively) for a fixed field mag- 
nitude of 700 mT. The data in (a) correspond to a rotation of 
the field around the /13 axis for the lowest-energy Zeeman-split 
energy-level doublet from an initial direction along fi2 to a fi- 
nal direction along hi, while plots (c) and (e) correspond to 
rotations from /13 to hi and /13 to /12, respectively. The color 
scale is identical to that used in Fig. 1. The g-factor of the 
first Zeeman-split doublet is shown in (b), (d), and (f) as 
a function of the rotation angles and is compared with the 
ellipsoidal dependence expected from theory. 

quantum states in the same sample, both in the magni- 
tude of the principal-axis g-factors and in the direction 
of maximum field sensitivity. 

Random matrix theory makes quantitative predictions 
about the statistical distributions of the g-factors along 
each of the principal axes, and how these distributions 
should depend on the strength of the spin-oijJpit coupling 
parameter (A 2 ocl/r^ ). Within the theory of0, the distri- 
butions of #1, #2? and g 3 may be generated for any A, and 
the value of A corresponding to a particular nanoparti- 
cle is determined by matching to the measured quantity 
(91)^(92)^(93)' Here () represents an average over en- 
ergy eigenstates in a particle. We find A=1.8±0.1 for 



TABLE I. Principal- axis g-factors and principal- axis direc- 
tions for individual energy levels in samples Cu#l and Cu#2. 
<7i, #2, and #3 are the principal-axis g-factors. 0\ and <j>\ are 
the spherical polar coordinates for the hi principal- axis di- 
rection while 63 and 03 give the /13 direction. is measured 
with respect to the z direction (perpendicular to the nitride 
membrane), while <j> is measured with respect to the x-axis. 
All angles are in radians, hi can be determined as hi x /13. 
The g-factors of states 6 and 7 in Cu#l and 1, 6, and 7 in 
Cu#2 could not be measured accurately due to nonlineari- 
ties associated with avoided crossings between closely-spaced 
states. 



level 


9i 


92 


93 


61 


0i 


63 


03 



Cu#l/1 0.86±0.08 0.51=b0.05 0.38=b0.05 0.71 4.64 0.99 2.22 

Cu#l/2 1.51±0.05 1.06=1=0.04 0.45±0.03 0.94 4.30 0.77 0.43 

Cu#l/3 1.65±0.09 0.79=b0.08 0.30=b0.06 0.13 3.13 1.54 4.76 

Cu#l/4 1.5±0.1 1.0±0.1 0.3±0.1 1.56 0.39 0.51 1.98 

Cu#l/5 1.34±0.06 0.78=t0.07 0.71=b0.07 1.56 1.31 1.54 6.06 

Cu#l/8 1.1±0.2 0.6±0.2 0.3±0.2 0.20 3.20 1.44 0.97 



Cu#2/2 1.3=1=0.2 1.2±0.2 0.9±0.1 

Cu#2/3 1.8±0.2 1.4±0.2 1.2±0.2 

Cu#2/4 1.6±0.3 0.9=1=0.2 0.6±0.1 

Cu#2/5 1.9±0.1 1.0=1=0.2 0.5±0.1 

Cu#2/8 1.3±0.2 1.2±0.1 0.9±0.1 

Cu#2/9 1.9±0.1 1.4±0.1 1.2=1=0.2 



0.46 5.72 1.24 3.37 
1.49 5.85 0.08 2.33 
1.51 1.89 0.98 0.28 
1.39 5.02 0.20 2.42 
0.74 1.21 0.83 4.52 
0.65 3.96 0.92 0.78 



Cu#l and A=l.l=t0.1 for Cu#2. Because spin-orbit 
coupling is thought to be dominated by scattering from 
defects and surfaces, it is not surprising that we find 
slightly different values of A for two particles made from 
the same materiallij. Given these A, the average val- 
ues for the g-factors along the principal axes, calculated 
assuming that the contributions to Zeeman splitting 
from orbital moments can be neglected relative to spin 
momentsQ'El are predicted to be (gi)=1.25, (#2) =0.76, 
and (g 3 )= 0.52 for Cu#l, and for Cu#2 ( gi )= 1.59, 
(#2)= 1.12, and (#3)= 0.96. The experimental aver- 
ages calculated from Table 1 are (gi)= 1.3=t0.3, (#2) = 
0.8±0.2, and (g 3 )= 0.4±0.2 for Cu#l and {g 1 }= 1.6±0.3, 
(g 2 )= 1.2±0.2, and (g 3 )= 0.9±0.3 for Cu#2. Therefore, 
using one fitting parameter (A) per sample, our results 
are in excellent agreement with the predictions. 

The distribution of principal-axis directions can pro- 
vide insight as to whether the axes are preferentially 
aligned due to crystalline axes or sample shape, or 
randomly oriented, as would be expected for spins 
coupling to highly-oscillatory electron-in-a-box orbital 
wavefunctionsQ. We tested whether the principal axes 
that we measure show a tendency toward any preferred 
direction by plotting the axes on polar plots (Fig. 3). 
The angles and <p are spherical polar coordinates, with 
measured relative to z, the unit vector perpendicular 
to the plane of the sample's silicon nitride membrane. 
The circles of constant are spaced so that equal areas 
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FIG. 3. Distributions of principal-axis directions, (a), 
(b), (c). Data points represent the hi, fi2, and fi3 princi- 
pal-axis directions, respectively, with circles for sample Cu#l 
and squares for Cu#2. The numbers identify the energy levels 
(see Table I) , in order of increasing excitation energy. and 4> 
are spherical polar coordinates, with measured relative to z, 
the direction perpendicular to the plane of the silicon nitride 
membrane. The radii of the circles of constant are scaled so 
that the areas inside the circles are proportional to the solid 
angle of the sphere with polar angle less than 0. 

of solid angle on the hemisphere are mapped to equal ar- 
eas in Fig. 3. The principal axes for different quantum 
states clearly do vary widely, and there does not appear 
to be any clustering in particular directions. 

In summary, we have made the first measurements of 
the angular dependence of the ^-factors for spin Zeeman 
splitting for individual quantum states in quantum dots. 
The results verify the predictions that spin-orbit coupling 
can cause the spin properties of quantum states to de- 
pend sensitively on the direction of an applied magnetic 
field. When we compare the ^-factors between different 
quantum states in the same nanoparticle, there are large 
variations both in angular dependence as well as magni- 
tude, supporting the picture that the strong anisotropics 
arise from the intrinsic random variations of quantum- 
mechanical wavefunctions. 
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